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ABSTRACT 

Brightest cluster galaxies (BCGs) in the cores of galaxy clusters have distinctly different 
properties from other low redshift massive ellipticals. The majority of the BCGs in cool-core 
clusters show signs of active star formation. We present observations of NGC 4696, the BCG 
of the Centaurus galaxy cluster, at far-infrared (FIR) wavelengths with the Herschel space 
telescope. Using the PACS spectrometer, we detect the two strongest coolants of the interstel- 
lar medium, [CiiJ at 157.74 /im and [Oi] at 63.18 //m, and in addition [Nii] at 121.90 /zm. The 
[Cn] enaission is extended over a region of 7 kpc with a similar spatial morphology and kine- 
matics to the optical Ha emission. This has the profound implication that the optical hydrogen 
recombination line. Ha, the optical forbidden lines, [Nn] A 6583A, the soft X-ray filaments 
and the far-infrared [Cn] Une all have the same energy source. 

We also detect dust emission using the PACS and SPIRE photometers at all six wave- 
bands. We perform a detailed spectral energy distribution fitting using a two-component mod- 
ified black-body function and find a cold 19 K dust component with mass 1.6x10^ Mq and a 
warm 46 K dust component with mass 4.0x10^ Mq. The total FIR luminosity between 8 fim 
and 1000 /zm is 7.5x10^ Lq, which using Kennicutt relation yields a low star formation rate of 
0.13 Mq yr~^ This value is consistent with values derived from other tracers, such as ultravi- 
olet emission. Combining the spectroscopic and photometric results together with optical Ha, 
we model emitting clouds consisting of photodissociation regions (PDRs) adjacent to ionized 
regions. We show that in addition to old and young stellar populations, there is another source 
of energy, such as cosmic rays, shocks or reconnection diffusion, required to excite the Ha 
and [Cii] filaments. 



1 INTRODUCTION 

Clusters of galaxies offer us a unique opportunity to study astro- 
physical components on widely differing scales. Intensive theo- 
retical and observational efforts have revealed that these compo- 



nents are closely tied to one another. The megaparsec-scale intra- 
cluster medium (ICM) is a hot plasma emitting bremsstrahlung X- 
ray radiation. The central regions of this plasma (< 300 kpc) have 
high electron densities; hence the numerous observed clusters with 
peaked surface-brightness profiles. Such high gas densities imply 
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rapid gas cooling (<c I/Hq) leading to a steady, mass deposition 
rate (tens to a few hundreds Mq yr"') (e.g. Fabian|l994 1 and itera- 
tively increasing central densities. 

However, this picture is far from complete. High-resolution 
XMM-Newton spectral and Chandra imaging of cool-core (CC) 
clusters (e.g., Tamura et al.'200r;'Peterson et al.'2001"Kaastra et al.' 
[200T] [Peterson et a l. 2003, Xu et al. 2002; Sakelliou et al. 2002 
[Sanders et al.|20 08 Hudson et al. 2010) have shown very little cool 
gas in their cores with the central temperatures no lower than a fifth 
of the virial temperature. Heating by clustercentric active galactic 
nuclei (AGN) through outflows (possibly in conjunction with other 
processes like heat conduction, cosmic-ray heating and convection) 
appears to be an irreplaceable heating mechanism that prevents the 



gas from undergoing excessive cooling (e.g. Binney & Tabor|1995| 



Churazov et al.|2002[[Roychowdhury et al.|2004[[Voit & Donahue[ 



20051). Several studies have established morphological, statistical 



and physical correlations between X-ray properties of CC clusters 
and AGN at their centers (Burns 1990, BTrzan et al. 2004; Rafi'erty 
[et al. [20061 [Peterson & Fabian 2006 ; Mittal et a l.|2009„McNamara 



[& Nulsen[[2007[ [Dunn & Fabi an 2006, EdwardTeFaTIgOOT) that 
support active galactic nuclei as the primary heating candidates. 

Despite AGN and other plausible sources of heating which 
stop the ICM from cooling catastrophically, recent observations of 
many (~ 40) brightest cluster galaxies (BCGs) show that the in- 
tracluster medium gas is cooling but at a much suppressed level. 
They appear to have substantial amount of cold gas, of which a 
small fraction is forming stars (e.g., [Johnstone et al.|1987[[McNa-[ 
Edge & Fr ayer|2003| [O'Dea et al.[[2004[|2010| 



mara et al.[[1989 



2008i[Hicks & Mushotzky|2005|[Mittaz et al.|2001|[Allen|1995^ . It 



seems that the cold gas is able to somehow survive for long periods 
of time. 

A significant number of cool-core clusters show optical line- 
filaments, such as in A 426 (Perseus, [Co"nseUce et al.[200lT ), A 1795 
( [Crawford et al.|2005^ . [Crawford et al.| ( |1999^ did an optical study 
of > 250 dominant galaxies in X-ray selected galaxy clusters and 
showed that about 25 % of them have Ho- emission-line in their 
spectra with intensity ratios typical of cooling flow nebulae (also 
see [Hu et al.] [l985 , Heckman et al. 1989[ ). Amidst these exciting 
discoveries of Ho-, Canning et al.[ ( [2010^ have detected the optical 
coronal line emission [FeX] X 6374A in NGC 4696 using VIMOS 
spectra, implying yet another comp onent of gas at temperatures in 
the range (1 to 5) x 10* K (also see Oegerle et al.j20 0p,'Bregman' 
eraL[|2006) . Warm molecular Hj at ~ 2000 K ^affe & Bremer 
T997[ [Donahue et al. 2000, Edge et al. 2002; Hatch et al. 2005), 



cold molecular H2 at (300-400) K (Johnstone et al..2007j and cold 
CO gas at few tens of kelvin pdge|2001[ [Salome et al.|2008| l also 
exists in the cores and also regions overlapping with the extended 
optical filaments. Edge (2001 ) detected 16 CC clusters with IRAM 
and JCMT, implying a substantial mass (lO'""-') Mq of molecular 
gas within 50 kpc radius of the BCGs. Similarly, [Salome & Combes[ 
([2OO3J studied 32 BCGs in cool-core clusters wifli the IRAM 30 m 
telescope and found gas masses between 3 x 10^ Mq to 4x 10'° Mq. 

From the above observations of CC clusters it is clear that 
there are different components of gas at almost every tempera- 
ture in the range varying from the virial temperature of the ICM 
(lO' K) to the temperature of the star-forming molecular clouds 
(10 K). An important result from the point of view of energetics of 
the observed warm molecular hydrogen emission and also the opti- 
cal coronal emission line is that these emissions are very likely due 
to the gas being reheated rather than cooling out of the ICM (I Don-[ 
[ahue et aT . 2000; Ca nning et al.|2010) . The molecular Ht emission, 
for example, is too bright to arise from gas cooling out of the ICM. 



Based on the measured H2 line luminosity and assuming the frac- 
tion of cooling in the line emission, t; = 2 % - 10 %, the estimated 
mass cooling rate is two orders of magnitude larger than the mass 
deposition rate predicted from the X-ray surface brightness profile. 
Understanding the details of how the mass and energy transfer oc- 
curs is crucial. To that end a key component is the molecular gas 
and dust at < 60 K, whose natural emission is accessible by Her- 
schel. 

In this paper, we study NGC 4696, the brightest galaxy of the 
Centaurus cluster of galaxies (Abell 3526). This work is a part of 
a Herschel Open Time Key Project devoted to study cold gas and 
dust in 11 BCGs. Preliminary results of this study can be found 
in [Edge et al. (2010b ) andlEdje et al. (2010a). NGC 4696 is at a 
redshift of 0.01016 ( ^Postman & Lauer„1995^ corresponding to a 
radial velocity of 3045 km s"'. Owing to its proximity, this cluster 
has been a subject of numerous studies (e.g. 'Fabian et alj 1982[ 
Lucey et al.|1986(|de Jong et al. 1990, Sparks et al. 1989; I O'Dea 
et al. 1994, Allen & Fabian 1994, Sparks et al. 1997; Laine et al. 
2003, ^Crawford et al.n2005( ,Taylor et al.,|2006( ,Johnstone et al 



2007[[Canning et al.|2010[[Farage et al.|2010^ . The Chandra X-ray 



observations show a bright dense core with plume-like structures 
that spiral off clockwise to the northeast. It is a classic cool-core 
cluster based on its short central [at 5 kpc (0.004 i^50(Q] cooling 
time (< 0.5 Gyr). However, the expected gas mass deposition rate 
within the cooling region is relatively small, ~ a few tens of solar 
masses per year ( Fabian et al. 1 982 ; H udson et al.[2010^[Ikebe et al.[ 
[l999[ [Sanders & Fabian 2002 ). Detailed XMM-Newton Reflection 
Grating Spectrometer (RGS) observations indicate the presence of 
cool gas at the center of the cluster with temperature in the range of 
(0.35 to 3.7) keV. Optical observations reveal bright line-emitting 
filaments in [Nii] A 6583A and Ha, which were first discovered by 
,Fabian et al.| l |1982) and have been mapped more recently using the 
EMMI instrument on the 3.58 m New Technology Telescope with a 
much higher resolution by jCrawford et al.] ( [2005[ ). These filaments 
extend out in a similar manner to the spiral structure seen in the X- 
ray. They also show a remarkable spatial correlation with the dust 
features, in particular the dust lane seen looping around the core of 
the NGC 4696 ( [Sparks et al.[T989[[Laine et al.|2003^ . 

NGC 4696 is a host to the low-power steep-spectrum FR I 
radio source, PKS 1246-410 (e.g. [Taylor "eTaT, 2006). The total 
radio luminosity of the radio source is 9.43 X 10*" erg s"' {Mil- 



[tal et al.[[2009j ). The radio core coincides well with the duo-core 
optical nucleus revealed by the high-resolution HST imaging. The 
radio emission, which shows a one-sided jet oriented to the south 
on small (< 30pc) scales that into lobes oriented east-west on kpc- 
scales l |Taylor et al.|2006l ). The total extent of the lobes along the 
east- west direction is about 10 kpc after which both the lobes even- 
tually turn south. The radio plasma is clearly interacting with the 
hot X-ray gas and seems to be correlated with the X-ray cavities. 

We describe the Herschel observations and data reduction in 
Sect. |2] and the results in Sect. [3] We present an analysis of the 
AGN contribution and kinematics in Sect. [4| and a detailed mod- 
eling of the photodissociation region in Sect. [5] We give our final 
conclusions in Sect. [6] and a summary of this study in Sect.|7] We 
assume throughout this paper the ACDM concordance Universe, 
with Ho = 71 hji km s"' Mpc"', = 0.27 and CIa = 0.73. This 
translates into a physical scale of 1" = 0.2 kpc at the redshift of 
NGC 4696. 



' R500 is the radius within which the average mass density of the cluster is 
500 times the critical density of the Universe. 
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2 HERSCHEL OBSERVATIONS AND DATA ANALYSIS 



We used the ESA Herschel space observatory ^Pilbratt et al. 



using PACS ( |Poglitsch et al.|2010| l and SPIRE (Grif fin et al. 



2010 



2010 



to study the cold phase (< 60 K) of the interstellar and intraclus- 
ter medium. The aim of these observations was to understand the 
details of mass and energy transfer between the different phases 
of gas. In particular, we study photodissociation regions (PDRs) 
which are relatively small in volume-filling factor but owing to high 
densities (and thermal pressures) as compared to that of the average 
interstellar medium (ISM) dominate the radiation of a galaxy (e.g. 
IHoUenbach & Tielens|1999t . PDRs define the emission character- 
istics of the ISM and star formation regions in a galaxy. They pri- 
marily comprise molecular hydrogen, ionized carbon, neutral oxy- 
gen, CO and dust, such as silicates, silicate carbides, poly aromatic 
hydrocarbons (PAHs) etc. Detailed PDR descriptions and chemical 
processes can be found in (e.g.|Hollenbach & Tielens|1999[|KauFl 



[man et al.|1999i|R611ig et al.|2007 1 



2.1 PACS Spectrometry 

We observed the two primary coolants of the neutral ISM, the 
[Cn] line at 157.74 //m and the [Oi] line at 63.18 yum, along with 
[Oib] at 145.52 yum, [Sii] at 68.470 yum, [Nii] at 121.90 ^im. and 
[Om] at 88.36 /jm with the PACS spectrometer on Herschel. The 
^P3/2 — >^Pi/2 fine-structure emission line of [Cii] is very often the 
brightest emission line in galaxy spectra, followed by the [Nii] lines 
at 122 yum ('Pj -»''Pi) and 205 /jm (^Pi ->^Po). Ionized nitrogen 
and oxygen are mainly produced in warm ionized medium, such as 
Hn regions, and so the [Nii] 122 yum and the [Oiii] 88.36 fim lines 
can be used to obtain constraints on the fractions of ionized and 
neutral media. 

All lines were observed in the line-spectroscopy mode with 
the chopping and nodding implementation to subtract the rapidly 
varying telescope background and dark current. The details of the 
observational parameters are summarized in table [T] We used the 
large chopper throw throughout, such that the mirror chops two 
regions of sky 6' apart alternately. While [Oi] and [Cn] observa- 
tions were made in the raster-mapping mode, all the other line ob- 
servations were made in the pointed mode with a field of view of 
47" x47". The PACS spectrometer contains an image slicer unit for 
integral field spectroscopy. The slicer transforms a 5x5 pixel focal 
plane image into a linear array of 25 spatial pixels termed spaxeh, 
each 9.4" x 9.4" in size. This technique ensures high sensitivity 
to detecting weak emission lines and provides simultaneous spec- 
tra of all extended emission in the field of view. The signal from 
each spaxel then goes through the grating assembly resulting in 16 
spectral elements. The PACS spectrometer has two channels, red 
(51-105) yum and blue (102-220) yum, which can be operated at the 
same time. This allowed us to perform sequential observations of 
both [Oi] and [Cii]. The raster mapping consisted of a 3 by 3 array 
of overlapping single pointings (3 raster lines and 3 points per line) 
separated by 23.5" in both the directions, along and orthogonal to 
the raster lines. 

We applied the standard pipeline routines described in the 
PACS data reduction guideline (PDRG) to process the data from 
their raw state (level 0) to a fully-calibrated state (level 2), us- 
ing the Herschel Interactive Processing Environment (HIPE) (Ott 
|2010[ l version 3.0 CIB 1475. The data-reduction steps included re- 
moving cosmic -ray glitches and correcting for the intrinsic non- 
linearities. We used the ground-based flat-field estimates to cor- 
rect for the pixel responses, which are known to overestimate the 



fluxes. We followed the PACS spectroscopy performance and cal- 
ibration (PSPC) memo to correct the fluxes by applying factors of 
1.3 for the blue band and 1.1 for the red band. The penultimate 
final PACS data product is in the form of two 5 x 5 x n cubes, cor- 
responding to the two nodding positions, where n is the number 
of wavelengths present. This product is further rebinned in wave- 
length in accordance with the actual wavelengths present and the 
user- requested resolution (see Sect. |2.1.T) . The PACS spatial reso- 
lution varies from 1 1.3" for the [Cn] 157.74 yum line to 4.5" for the 
[Oi] 63.18 //m line. 



2.1.1 Over- and up-sampling factors 

As mentioned above, the 5 x 5 x n cubes are further rebinned in 
wavelength by applying two factors called 'oversample' and 'up- 
sample'. The oversample factor increases the number of bins over 
which the data are averaged to create a spectrum and the upsam- 
ple factor specifies the shift in wavelength in the units of the bin- 
width. The default values, oversample=2 and upsample=l, corre- 
spond to the signal sampled at the Nyquist-Shannon rate, yielding 
the instrument spectral resolution at that wavelength. Increasing the 
upsample factor, for example, to upsample=2 yields two spectra, 
both adhering to the instrument resolution but shifted by half a bin- 
width with respect to each other An upsample factor greater than 
unity may be used to identify features narrower than the instru- 
ment resolution in the line profile but at the expense of increased 
noise. Lastly, the combination oversample=l and upsample=l un- 
dersamples the data to render a spectral resolution a factor of two 
coarser than the original one, hence affecting the full-width half 
maximum (FWHM) of the fitted line. We investigated various com- 
binations of oversample and upsample factors for each of the lines 
to optimize the signal to noise ratio and the spectral resolution. 



2.1.2 Line Flux Estimation 

In the case of emission well-centered on the central spaxel, a 
wavelength-dependent point-source correction factor has to be ap- 
plied to the measured line flux-density to account for the fraction of 
the beam that is lost to the neighbouring spaxels (beam spillover). 
An alternate way of recovering the flux in the entire beam is to co- 
add the spectra obtained in different spaxels but this method suffers 
from the risk of line distortion. A yet alternate method is to keep 
the spectra separate and add only the fluxes contained in the spaxels 
with significant point source contribution. This method works well 
also for slightly extended sources provided that the emission is fully 
contained within the central few spaxels of the pointing. However, 
in case of extended emission spread over more than a single point- 
ing (47" x47"), we used the HIPE inbuilt task 'specProject' for ob- 
taining the integrated line flux. This algorithm involves a projection 
of all the positions in the input cube onto the sky. It calculates a reg- 
ular RA/Dec grid with a default spaxel size of 3" and projects the 
raster positions onto it. The details of this algorithm can be found 
in the spectroscopy pipeline described in the PDRG documenta- 
tion. Thereafter, we conducted 'aperture photometry' by placing 
a rectangular box around the visible emission and integrating the 
flux to obtain a spectrum. Since none of our observations were per- 
formed in dithered mode, the main necessity of using 'specProject' 
was to obtain a projection of the multiple single pointings on the 
sky plane. Since the individual pointings overlap by a considerable 
amount, simply adding the fluxes in all spaxels contributing above 
a certain threshold results in an overestimation of the line emission. 
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Table 1. Herschel PACS spectroscopy observational log of NGC 4696 at a redshift of 0.01016. All the lines were observed in the line spectroscopy mode and 
on the same day: 30"" Dec. 2009. 



Line 


Peak Rest A 


ObsID 


Duration 


Bandwidth 


Spectral FWHM 


Spatial FWHM 


Mode 




(fim) 




(s) 


(//m) 


(kms-') 


(/jm) 


(km s-') 






Ol 


63.180 


1342188700 


6912 


0.266 


1250 


0.017 


79 


4.6" 


3x3 raster, step size 23.5" 


CII 


157.74 


1342188700 


3096 


1.499 


2820 


0.126 


237 


11.5" 


3x3 raster, step size 23.5" 


Nil 


121.90 


1342188701 


3440 


1.717 


4180 


0.116 


280 


8.9" 


pointed 


Olll 


88.36 


1342188701 


3680 


0.495 


1660 


0.033 


110 


6.4" 


pointed 


Olb 


145.52 


1342188702 


3440 


1.576 


3215 


0.123 


250 


10.5" 


pointed 


Sil 


68.47 


1342188702 


3840 


0.218 


945 


0.014 


62 


4.9" 


pointed 



We used the Levenberg-Marquardt minimization routine (fLev-| 
|enberg| 1944[ |Marqua rdt 1963) to fit a model to an observed spec- 
trum. The model comprised a gaussian for the line emission and a 
polynomial of order (or 1) for the pseudo-continuunj^ baseline. 
A positive line detection was based in terms of the signal-to-noise 
ratio (SNR). The SNR was defined as the ratio of the line peak to 
the standard deviation of the data about the fitted model. 

2.2 PACS Photometry 

The PACS photometric observations were taken in large scan 
mapping mode in all three bands of the photometer, blue- 
short (BS) (70 fim), blue-long (BL) (100 t^m) and red (R) (160 ^m) 
using the medium scan speed (20"s"'). The scan maps comprised 
18 scan line legs of 4' length and cross-scan step of 15". Each ob- 
servation had a "scan" and an orthogonal "cross-scan" direction 
and we calibrated the corresponding data separately before com- 
bining them into a single map of 9' x 9'. The resulting maps have 
a resolution of 5.2", 7.7" and 12" at 70 yum, 100 yum and 160 lum. 
The PACS photometer performs dual-band imaging such that the 
BS and BL bands each have simultaneous observations in the R 
band so we have two sets of scans in the R band. 

For PACS photometry, we used the HIPE version 
3.0 CIB 1475 to reduce the data and adopted the PDRG to 
process the raw level-0 data to calibrated level 2 products. We 
employed the official script for PACS ScanMapping mode with 
particular attention to the high pass filtering to remove "V-y'/" 
noise. We used the 'HighPassFilter' method to remove the large 
scale (> 80") artefacts. The target BCG and other bright sources in 
the field were masked prior to applying the filter. The size of the 
mask was chosen to be less than the filter size so as to minimize 
any left-over low-frequency artefacts under the masks. We used a 
filter size of 20 readouts for the BS and BL bands and 30 readouts 
for the R band and a mask radius of 20" for the BS and BL bands 
and 30" for the R band. Finally the task 'photProject', was used 
to project the calibrated data onto a map on the sky in units of 
Jy pixeL'. The "scan" and "cross-scan" maps were then averaged 
together to produce the final coadded map. 

2.3 SPIRE Photometry 

The SPIRE photometry was performed in the large scan map mode 
with cross-linked scans in two orthogonal scan directions. The pho- 
tometer has a field of view of 4' x 8', which is observed simulta- 
neously in three spectral bands, PSW (250 //m), PMW (350 fim) 

^ Presently, there are offsets in the continuum levels between the different 
spaxels preventing the extraction of reliable continuum data. 



and PLW (500 yum) with a resolution of about 18", 25" and 36", 
respectively. 

For SPIRE photometry, we used a newer HIPE version 
4.0 CIB 1432 and the standard HIPE pipeline for the LargeScan- 
Map observing mode and the naive map-maker. The pre-processed 
raw telemetry data were first subject to engineering conversion 
wherein the raw timeline data were converted to meaningful units, 
the SPIRE pointing product was created, deglitching and temper- 
ature drift correction were performed, and maps were created, the 
units of which were Jy beam"' . 



3 RESULTS 

3.1 Line Detections 

Of the lines observed, we detected [Oi], [Cii] and [Nn] in 
NGC 4696 (Fig.[T). We did not detect [Sii], [Oni] and [Om] emis- 
sion lines at a 3cr level of 3.0 x 10 W/m^ for [Sii] and [Om], and 

1.2 X 10 W/m^ for [Oib]. The line parameters are summarized in 
tableU 

The [Nil] emission line at 121.9 fim has been detected only 
in the central 9.4" spaxel. The [Oi] emission line spectrum at 
63.18 yum, although noisy, has clearly been detected. Shown in the 
top left panel of Fig.[T]is the continuum-subtracted spectrum inte- 
grated over a central region of ~ 15" in diameter (3 kpc). Hence, the 
[Oi] emission is slightly extended. The spectrum is obtained with 
oversample=2 and upsample=l followed by the projection tech- 
nique, 'specProject', described in Sect. |2.1.2| 

The [Cii] continuum-subtracted integrated line emission at 
157.74 fim is shown in the top right panel of Fig. [T| It is obtained 
using oversample=l and upsample=4 and projecting the two nods 
onto the sky. Though oversample= 1 results in poorer instrumental 
resolution, it helps in beating down the noise and upsample=2 or 

4 helps recover some of the lost information. The two-dimensional 
distribution of [Cii] emission is shown in the multiwavelength col- 
lage in Fig. |2] which contains an optical Ho- image from NTT, 
an X-ray image from Chandra and 1.4 GHz radio contours from 
VLA/NRACp] The Chandra datasets used to create the upper left 
panel of Fig.[T]were Obslds 504, 4954, 4955 and 5310, which were 
first presented in |Sanders & Fabian] ( |2002^ and |Fabian et al.|f2002| 
[2005 ). The archival datasets were reprocessed with ciao 4.3 and 
CALDB 4.4.0 to apply the appropriate gain and charge transfer ineffi- 
ciency correction, then filtered to remove photons detected with bad 
grades and cleaned of ffare periods. The final cleaned exposure time 

' The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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T 





Velocity [km/s] (OI line) 



A 



1000 

p 



111 



63.60 63.65 63.70 63.75 63.80 63.85 63.90 63.95 64.00 64.05 

Wavelength [//m] 



Velocity [km/sl (CII line) 

-2000 -1500 -1000 -500 500 1000 1500 




158.2 1584 1.58.6 1.58.8 159.0 159.2 1.594 1.59.6 159.8 160.0 160.2 160.4 

Wavelength \ji m] 



Velocity [km/s] (NIT line) 

-2500 -2000 -1.500 -1000 -.500 .500 1000 1500 2000 2500 




122.0 122.5 123.0 123.5 124.0 

Wavelength [/im\ 

Figure 1. The forbidden far-infrared line detections in the center of NGC 4696 made with the Herschel PACS instrument. The lines are spatially integrated: 
the [Oi] 63.18 fim line (upper left panel, marginally extended), the [Cii] 157.74 fim line (upper right panel, extended) and the [Nn] 121.9 /im (lower panel, 
point-like). 

Table 2. Estimated parameters for the far-infrared forbidden lines in NGC 4696. Also given are the upper limits on the line flux for the non-detections. The 
total extent of the [Cii] and [Oi] emission was estimated based on visual inspection. 



Line 


A dim) 


Offset (km s"') 


FWHM (kms-') 


Line Flux 


Total Extent 






^bcg 




Obs. 


Intrinsic 


(10"'**) W/m2 




[Oi] 


63.827±0.005 


-15±28 


-52±28 


228±54 


218±35 


57.6±7.7 


15" 


[Cii] 


159.370±0.003 


18±11 


-22± 1 1 


410±13 


335±16 


174.7±3.1 


35" 


[Nil] 


123.171±0.009 


48±40 


8±40 


452±54 


351 ±69 


24.9±1.7 


< 9.4" 


[Oiii] 












< 3 




[OlB] 












< 2 




[Sii] 












<4 





was 20L6 ks. Each of the final cleaned events files was then repro- 
jected to match the position of the 4954 observation. A combined 
image was produced by summing images in the (0.5 to 1.0) keV 
band extracted from the individual reprojected datasets. This image 
was then corrected for exposure variation by dividing the summed 
image by summed exposure maps created for each dataset. 

The total extent of the [Cn] emission is about 35" (7 kpc) in di- 
ameter. The [Cii] image exhibits a northwest-southeast elongation, 
which in conjunction with the optical and X-ray emission maps of 
the spiral filaments, is clearly the overall elongation direction asso- 
ciated with NGC 4696. The PACS continuum dust emission shows 
similar asymmetry about the center of the galaxy (Fig.[3|. The bulk 



of the [Cii] emission appears to be displaced to the west relative 
to the Hq' core. The centroid of the [Cii] emission is indeed offset 
to the west of the Ho' core by about 5" (half spaxel size). The re- 
ported lo" absolute pointing accuracy of Herschel is 2" for pointed 
observations and larger for scan map observations. As a check, we 
analysed the [Cn] line observations with the latest official version of 
HIPE (6.0.2055) with updated calibration files and found a similar 
distribution of the [Cii] emission with a similar offset from the Ho' 
core. While the offset between the [Cii] and Via emission may also 
be real and reflect inhomogeneities in the ISM properties or differ- 
ent excitation mechanisms, given the similarity between [Cii] and 
Ha it is more likely that the offset is due to the pointing uncertainty. 
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Figure 2. Upper panel: The Chandra X-ray surface-brightness map in the soft energy band (0.5 - 1) keV (left), the Gemini optical Hff emission (middle) 
[Crawford et"ar](2005) and the Herschel far-infrared [Cn] 157.74 /.im emission (right). ; Lower panel A multiwavelength RGB image of NGC 4696 with 
contours. Red: optical Ha line emission. Green: X-ray emission. Blue: far-infrared [Cii] emission. Black contours: VLA 1.4 GHz radio contours (NRAOA'LA 
image archive). Each side measures an arcminute corresponding to a linear scale of 12 kpc. 



3.2 Aperture Photometry 

We conducted aperture photometry using a variety of different 
methods with the aim of attaining robust flux measurements and re- 
liable uncertainties. These methods included routines inbuilt within 
HIPE, GAIA (Graphical Astronomy and Image Analysis Tool) etc. 
The spread in the flux estimates was within the random noise asso- 
ciated with the measurements. 

For PACS flux densities, small aperture corrections were ap- 
plied as outlined in the PACS ScanMap release note. Further flux 
calibration on the derived flux densities was performed to account 
for the known overestimation introduced by the ground-based flux 
calibration by factors of 1.05, 1.09 and 1.29 in the BS, BL and R 
bands respectively. The PACS absolute flux accuracy is within 10 % 
for BS and BL, and better than 20 % for R. The uncertainties are not 
believed to be correlated due to the BS and BL bands being taken at 
dififerent times and the R band using a different detector. Since we 
used a newer HIPE version to obtain SPIRE images, the known flux 
calibration offset was accounted for during the data conditioning 
and hence no multiplicative calibration factors were required. The 
SPIRE absolute flux accuracy is within 15 % for all three bands. 
The flux measurements froin Herschel and other instruments are 
given in table|3] 

Our Herschel flux estimate disagrees with the Spitzer MIPS 
160 ^lm measurement of Hoffer et al. ^2011^, th ough it is consistent 
with the measurement of |Kaneda et al.| ( |2005^ . A systematic error 
in the performance of either of the two instruments can be ruled 
out since while at 70 yum the Spitzer flux estimate is higher than the 



Herschel estimate, at 160 ^m the situation is reversed. In addition, 
we also ruled out the discrepancy being due to the difference in the 
spectral response functions of Herschel PACS and Spitzer MIPS 
photometers. We used the above fitted SED of NGC 4696 as the 
'true' source model and calculated the predicted flux density for 
each instrument as the weighted average, where the weights were 
determined from the spectral response function of each instrument. 
The Herschel PACS and Spitzer MIPS predicted flux densities at 
70 jim and 160 yum estimated this way are consistent with each 
other to better than 5 %. Since the Spitzer MIPS 160 pm flux shows 
a spread between the measurements from |Kaneda et al.] ( |2005[ l and 
|Hoffer et al.|pOI \\ , we give preference to the Herschel estimates, 
especially since Herschel is designed to operate at higher sensitivity 
at FIR wavelengths. 



3.3 Dust SED 

The dust emission froin NGC 4696 has been detected in all 
six PACS and SPIRE bands. In Fig. [5] we show the PACS im- 
ages at 70 /jm, 100 yum and 160 yum and the SPIRE images 
at 250 yum, 350 yum and 500 yum centered on NGC 4696. The 
PACS photometry data reduction pipeline described in Sect. |2.2| 
produces images with units Jy pixeL'. We converted the units 
from Jy pixeL' to Jy beam"' by using the conversion factor 
7r/(41n2)(0beam/Spix)^ pixel beam"', where Sbcam is 6", 7" and 12" 
and epi„ is 3.2", 3.2" and 6.4" for the 70 yum, 100 //m and 160 yum 
bands, respectively. 
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Figure 3. Photometry images of NGC 4696. Left: PACS colour images at 70 fim (blue), 100 /jm (green) and 160 yum (red) with a resolution of 6". 7" and 
12" respectively. Center. PACS colour images combined with the same smoothing gaussian of FWHM 12". Right: SPIRE colour images at 250 /im (blue), 
350 fim (green) and 500 /jm (red) with a resolution of 18", 24" and 38" respectively. The image units are Jy/beam and each side measures 7'. 
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Figure 4. The spectral energy distribution (SED) of NGC 4696. Left: The far-infrared SED obtained from a modified two-component blackbody fit to the 
Herschel and Spitzer data, designated as 'FIT' (filled black diamonds). Right: The complete SED extending to radio frequencies on one end and optical and 
X-ray frequencies on the other, collectively designated as 'MISC. 



The BCG emission seems to be slightly extended along the 
northwest-southeast direction. The origin of the infrared continuum 
emission is the interstellar dust grains in the host galaxy, which 
are heated by the absorption of photons from an ionizing source. 
Whether this ionizing source is starbursts (SB) or an AGN em- 
bedded in the host galaxy is usually a matter of debate but the 
AGN/SB diagnostics for NGC 4696 suggest the AGN contribution 
to be very low (Sect. pTTj ). We fitted the dust spectral energy distri- 
bution (SED) with a blackbody spectrum modified by a frequency- 
dependent emissivity. Following |Aravena et al.| ( f2008) , and refer- 
ences therein, the flux-density at a given frequency, v, is : 

Sy = [BATi) - i?,,(r,„b)] (1 - (1) 

(1 +z)' 

where By(T) is the Planck function at frequency v, and temperature 
T. B,,(rcmb) is the contribution from the cosmic microwave back- 
ground at Tcmb = 2.73(1 -I- z) K. n is the solid angle subtended by 
the source, n(f, where 6 is the source size in radians (here assumed 
to be the resolution of the PACS photometer at 70 jim). The last 
term in parenthesis is the modification to the blackbody radiation 
spectrum, where is the dust optical depth. 



dust absorption coefficient and /J is the dust emissivity index which 
we fixed to 2 . 



T,, = Ky 



DlSl 



(2) 



Here, Ky = 5.6 x (v/3000 GHzy m^ kg ' i Dunne et al. 2000 > is the 



Motivated by the findings of |Dunne & EalesV2001l), we fitted 
the infrared data with two modified blackbody functions, represent- 
ing a cold dust component and a warm dust component. This gives 
four model parameters: the temperature and mass of the cold dust 
component, Td c, Md c, and the temperature and mass of the warm 
dust component, T^ ,,, Mj „. We used the Levenberg-Marquardt 
minimization routine from Numerical Recipes to obtain the best 
fit model parameters. We conducted the fit using seven data points: 
the six PACS and SPIRE photometer data and the Spitzer MIPS 
data at 24 yum. The Spitzer MIPS flux density at 24 yum receives 
a contribution from the hot dust component as well as the passive 
stellar population. In order to use the 24 ytim data for the SED fit- 
ting, it is necessary to subtract the contribution from the passive 
stellar population. For this we fitted a powerlaw to the IRAC data, 
since the IRAC wavelengths are sensitive to emission from the old 
stellar population, and extrapolated the fit to 24 /im. The net flux- 
density after subtracting the contribution from the passive stellar 
population is 3.6 ± 3.0 mjy, which is significantly lower than the 
total observed value. It is this net flux-density that we assigned to 
the hot dust component for the SED fitting. 

The SED fitting yields a massive cold dust component with 
r,,., = (18.9 ± 0.7) K and Mj., = (1.6 ± 0.3) x 10*^ Mq and a 



8 Rupal Mittal et al. 



Table 3. Compilation of infrared flux-densities. 



A (pim) 


Instrument 


Aperture 


Flux (mjy) 


J.O 


TT? AC^ Qr\itvfiT-'^ 
IKrt^.. opiLZcl 


DJ 


JZo± 1 / 




iKi\^^ opiizer 


DJ 


1 o 1 -1- 1 n 

ly L± i\J 


J.O 


TT? AC Snif7pr^ 


DJ 


162+8 


8 


TR AC *snif7pr^ 


fi5" 
DJ 


1 nn+s 


12 


IRAS'^ 






Z4 


Mlr^o opuzei 


Dj 


Zj.O± J.U 


Z4 


Mir^i opitzer 




Z4±4 


25 


IRAS" 




< 28 


60 


IRAS*^ 




100±23 


70 


MIPS Spitzer^ 


35" 


165±33 


70 


MIPS Spitzer'' 


30" 


156±19 


70 


PACS Herschel'' 


36" 


114±11 


100 


PACS Herschel'' 


36" 


244±24 


100 


IRAS'^ 




830±148 


160 


PACS Herschel'' 


36" 


383±77 


160 


MIPS Spitzer^ 


40" 


256±51 


160 


MIPS Spitzer'' 


30" 


331±29 


250 


SPIRE Herschel'' 


38" 


208±30 


350 


SPIRE Herschel'' 


38" 


98±15 


500 


SPIRE Herschel'' 


38" 


38±6 



^ The Spitzer IRAC/MIPS data are from a sub-sample of the ACCEPT 
sample of galaxy clusters ( Hoffer et al.|201 1] in preparation). The error- 
bars reflect the absolute flux uncertainties. 



Kaneda et al. 
The IMS 



2005 
data 



2007 
T 



are trom NASA/IPAC Extragalactic Database 
jhttp://nedwww.ipac. caltech.edu i. 

" Herschel PACS/SPIRE data are from this work. The Herschel 
PACS/SPIRE errorbars reflect the absolute flux uncertainties: 10% at 
PACS BS and BL and 20% at R, and 15% at all SPIRE wavelengths. 



low-mass warmer dust component with T^^, = (46.0 ± 5.0) K and 
Md,^^ = (4.0 ± 2.8) X 10^ Mq. As also seen in other studies, the cold 
component dominates the mass budget. The far-infrared and the 
total SEDs, the latter extending from radio to X-ray frequencies, 
are shown in Fig.|4] We used the fit parameters to estimate the to- 
tal far-infrared luminosity of the dust between 8 /jm and 1000 yum, 
ipiR.tot. and found it to be (7.5±1.6)xlO^ Lq. For the analysis of 
the photodissociation regions described in Sect.|5] we used the FIR 
luminosity in the range 40 lum to 500 yum, which from here on is 
denoted by LpiR and is equal to (6.2±0.8)xlO* Lq. 

To assess the robustness of the derived parameters and the 
quantities derived from them, such as the FIR luminosity, we also 
fitted the above model using the average of the two MIPS and one 
PACS flux estimates at each 70 yum and 160 fim. The fit results 
lead to a cold dust component with Tdx = (18.2 ± 0.8) K and 
= (1-8 ± 0.3) X 10^ Mq and a warm dust component with 
= (43.0 ± 4.7) K and Md.„ = (9.1 ± 6.8) x 10' Mq. The 
inferred FIR luminosities are LFiR,toi = (7.9 ± 2.5) x 10* L© and 
LpiR = (6.6 ± 1.1) X 10** Lq. Hence, a -20% variation in the PACS 
flux densities results in about a 5% to 10% variation in the total 
cold dust mass and the FIR luminosities. This is well within the 
formal errorbars obtained from either of the two fitting methods. 



3.4 Gas-to-Dust Mass Ratio 



|0'Dea et al.| ( |1994} searched for CO in a sample of five radio galax- 
ies, including NGC 4696, using the Swedish-ESO submillimeter 
telescope. They reported nondetection in NGC 4696 and used the 
result to put limits on the gas mass. From these observations, the 



3cr upper limit on the mass in the molecular gas is estimated to 
be ~ 5 X 10** Mq. Using the dust mass calculated in Sect. 3.3 the 
(molecular) gas-to-dust mass ratio has an upper limit of ~ 325. 
Note that the upper limit on the total gas-to-dust mass ratio is likely 
larger since the total-to-molecular gas mass ratio is greater than 
unity and may be as high as 5 (or more) depending upon parame- 
ters such as the ionization radiation intensity and the total hydrogen 
column density. 

We also used the formalism laid out by |Wolfire et al.| ( [l990[ l, 
which uses the observed [Cii] emission to obtain the atomic gas 
mass, assuming that the [Cii] emission is optically thin: 



2.7 X 10' 



1 Mpc) 
10""' ergs s" 



^Cii 



A(Cii) 



10-" Wcm-2 

'^ /3xlO-^ 
^' 4Cii) 



(3) 



Mo 



where Z)l is the luminosity distance, Fq^ is the [Cii] line flux, 
A(Cii) is the cooling rate of ionized carbon and x(Cii) is the abun- 
dance of ionized carbon relative to hydrogen. Most of the [Cii] 
emission originates from the surface layers of the photodissocia- 
tion regions with the visual extinction. Ay < 4, depleting steadily 
with increasing Ay, and is the dominant of the three most important 
carbon-bearing species, C, Cn and CO ( HoUenbach et al.]|1991[ ). 
Hence the Cn abundance can be equated to the elemental abun- 
dance of carbon, which, for the simulations of the photodissocia- 
tion regions described in Sect. |5] equals 2.5 times the interstellar 
value, x(Cn) = 6.3 x 10"*. For gas temperatures higher than 100 
Kelvin (the excitation temperature of Cn is 92 K) and gas densi- 
ties between 100 cm"' and 1000 cm"', A(Cn) is in the range (10"^' 
- 10"^^) ergs s"' sr"' atom"'. Substituting these values into Eqn.|4] 
gives a gas mass in the range (0.5 to 5) xIO'Mq. This is consis- 
tent with the upper limit derived from the non-detection of CO by 
lO'Dea et al.l l pMt . 

The gas-to-dust mass ratio based on the above estimate of the 
gas mass is between 1 and 70. In comparison with the gas-to-dust 
mass ratios derived for other BCGs fEdge'2 00T||Edge et al.|2010al l, 
which have typical values of a few hundred, the ratio in NGC 4696 
is significantly lower. Note that the above method of determining 
the gas mass is not very robust for low n and GO, which is the case 
of NGC 4696 (see Sect.|5|, and the dust mass is very sensitive to the 
assumed dust temperature, so the derived range of the gas-to-dust 
mass ratio should only be treated as a crude estimate. 



3.5 Star Formation Rates 

Assuming the dust heating to be due to the young stellar population 
and high dust opacity in the star forming regions, the Kennicutt re- 
lation Kennicutt|I998^ can be used to convert the LpiR.toi (derived 
in Sect. \33\ into a star formation rate (SFR). Using the conver- 



sion, SFR (Mq yr"') = 4.5 x IO"'"'LFiR.toi (erg s"'), and the derived 
ipiR.tot! we obtained a star formation rate of ~ 0.13 M© yr"'. How- 
ever, we expect a non-negligible contribution also from the general 
stellar radiation field, including the old stellar population in the host 
galaxy, which will lower the estimated SFR. Hence, the above value 
is an upper limit. 

We also retrieved far ultraviolet (FUV) data froin the HST 
archive. The data (proposal ID 11681) were taken with the Ad- 
vanced Camera of Surveys Solar Blind Channel (ACS/SBC) pro- 
viding a field of view of 34.6" x 30.5". These observations were 
made with two long-pass filters, F150LP and F165LP, and we chose 
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S1.0 50.5 12:43:50.0 49.5 49.0 48.5 48.0 47.5 

Right ascension 

Figure 5. The ACS/SBC FUV F150LP observations of NGC 4696 with 
the HST showing a core surrounded by some diffuse low-brightness emis- 
sion. The extent of the FUV emission is very small relative to Ho- and [Cii] . 
Weak FUV emission is consistent with the low star formation rate derived 
in Sect. |3.5| and small values of the FUV intensity field. GO, described in 
Sect.H] 



the F150LP filter over F165LP to estimate the star formation rate 
(described below) due to its higher system throughput. 

We used a Starburst99 stellar library jLeitherer et al.||T999^ 
to model a spectrum of a young stellar population (whose old- 
est stars are 10' yr old) with a continuous star formation rate of 
0.1 Mq yr"'. We fed the output spectrum into the synphot syn- 
thetic photometry package to determine the expected flux in the 
F150LP band. The ratio of the extinction-corrected observed flux 
to the expected flux from synphot multiplied by the assumed SFR 
in Starburst99 gives a measure of the actual star formation rate. 
The FUV emission in the F150LP band, shown in Fig. [5] con- 
tains a bright core surrounded by low brightness emission. The 
low surface-brightness extended emission is most probably dom- 
inated by the passive stellar population. Hence, we derived an up- 
per limit to the SFR based on the assumption that the entire ob- 
served FUV emission arises due to a young stellar population and 
a lower limit based on the assumption that only the bright core 
is associated with an ongoing star formation. The measured flux 
using aperture photometry is 1.23 x 10"'^ ergs s"' cm"^ A"' if 
a region of radius ~ 1.5 kpc in considered, which encompasses 
the bright core as well as the diffuse emission surrounding it, and 
2.75 X 10"" ergs s"' cm"^ A"' if only the bright core of radius 
~ 0.5 kpc is considered. These values were corrected for Galac- 
tic extinction based on the measured E(B - V) = 0.113 from 
NED and the extinction law, Ry = 3.1[= A{V)/E(B - V)] ( |CardeUi| 
|et al.|1989] l. The extinction was calculated using the mean Ry de- 
pendent extinction law from jCardelli et aL| ( |1989[ ) in the UV and 
FUV regime at the pivot wavelength, A = 1612.236 A. The ex- 
tinction corrected fluxes are 2.80 x 10"'^ ergs s"' cm"^ A"' and 
6.27 X 10"" ergs s"' cm"- A"'. Note that there is likely to be inter- 
nal extinction within NGC 4696 and the SFR limits derived below 
may increase once this additional extinction is taken into account. 
Comparing the expected flux from the Starburst99 model to the ob- 
served corrected flux gives an upper limit of the FUV derived star 



formation rate of 0.08 Mq yr"' and a lower limit of 0.002 Mq yr"' . 
The upper value is in agreement with the upper limit derived from 
the far-infrared analysis. This is also consistent with the GALEX 
near-ultraviolet (NUV) observations indicative of a SFR consistent 
with zero and a 3<t upper limit of 0. 17 M© yr"' . The GALEX anal- 
ysis takes into account the UV-upturn from an evolved stellar pop- 
ulation. The UV excess is calculated by comparing the GALEX 
NUV emission to the 2MASi^ /f-band emission within an aper- 
ture of 35" and attributed to the ongoing star formation. 

The classical X-ray mass deposition rate at a radius of about 
35" from the X-ray peak in the Centaurus galaxy cluster (the pro- 
jected separ ation between the X-ray peak and the BCG is 6" , |Mittal| 



Sanders & Fabian 2008 



et al. 2009) is about (5 to 9) Mq yr"' (e.g. 

Hudson et al.||20f0l i. The cooling rates observed in the form of 
star formation rate derived from FIR and UV are smaller by about 
two orders of magnitude. Similarly, 'S anders & Fabian| ( [2008 l) used 
XMM-Newton RGS observations to study the X- ray emission from 
the core of the Centaurus galaxy cluster. They carried out detailed 
spectral fitting to measure the amount of gas cooling at the center 
and deduced an upper limit of 0.8 Mq yr"' below 0.4 keV. This 
is the standard "cooling-flow" discrepancy, which is resolved when 
we consider the AGN regulated feedback. The luminosity of the gas 
cooling in the center can be estimated using a cooling flow model 
in which the gas cools from the virial temperature to a minimum 
possible lower temperature (0.08 keV). The cooling luminosity is 
given by: 



5 Mo 
2 /imp 



-IrT ■ 



(4) 



where Mq is the classical mass deposition rate inferred from 
X-ray observations, njp is the mass of the proton, /j is the mean 
molecular weight (~ 0.65) and Tvir is the cluster virial temperature 
which for Centaurus is about 4 keV. Inserting the estimated M© 
and Tvi, gives a cooling luminosity in the range (4.5 to 8.0) lO''^ 
erg s"'. On the heating front, Rafferty et al. j2006| found the me- 
chanical power associated with the AGN -excavated cavities to be 
7.4x 10*^ erg s"' . Merloni & Heinz 1 2007^ estimated a similar value 
for the mechanical power ~ 7.8 x 10*^ erg s"'. Hence there is 
enough energy in the radio source to balance cooling of the hot 
gas and maintain equilibrium. 



4 ANALYSIS 

4.1 AGN contamination 

The interstellar dust is believed to be heated by two principle mech- 
anisms, by young stars and AGN. Hence it may be that some of the 
infrared emission detected is related to the AGN activity. In order 
to gauge the AGN contribution, diagnostics such as those used by 
[O'Dea et al.|p008) ; |Quirien et al.|(2008l i, such as the Spitzer IRAC 
4.5//m/3.6//m colour and [O m] /15007/H/3, may be used to discrim- 
inate between the presence of an AGN and star formation. 

We compiled the optical line emission ratio, [O m] ASOOT/H/}, 
for NGC 4696 available from |Lewis et^p003l l (< 0.4) and |Farage| 
|et al.|pOTO^ (~ 1). The two estimates are not necessarily contradic- 
tory since these are calculated over different apertures. Neither of 
the ratios indicates a presence of a dominant AGN contribution to 
the far-infrared emission. In addition to the low [O ni] ASOOT/Hp, 
the IRAC 4.5fim/3.6fim colour ratio (= 0.58) is consistent with a 



* http://irsa.ipac.caltech.edu/Missions/2mass.html 
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Figure 6. Gas kinematics. The measured line-of-siglit velocities (left panels) and line FWHM (right panels) in km s"' of the far-infrared [Cii] gas with a scale 
3" per pixel (upper panel) and the optical ^[^^-[Nii] A 6583A gas with scale 0.67" per pixel (lower panel). The velocity-dispersion maps have been corrected 
for the instrumental resolution (235 km s"' for [Cii] and 42 km s"' for Ha). The FIR [Cn] emission traces the optical Ha emission morphologically and 
kinematically. 



passive stellar population. The absence of the IR [Oiii] 88 jim spec- 
tral feature is another indication of negligible AGN photoionization 
( |SpinogUo|2()09f . 



4.2 Gas Kinematics 

Observations of NGC 4696 reveal interesting velocity structures 
and several studies have used data from different wavebands to 
study them (e.g. Sparks et al^lflggg] [TWTj |de Jong et aT1[T990ll. 
The most recent independent investigations by |Farage et al.H2010^ 
and |Canning et al.|j2011f [in prep.] use observations of Ha at sub- 
arcsecond resolution to map the kinematics of the gas. 

We used the projected map made using 'specProject' 
(Sect. |2.1.2| l and obtained a line fit for each individual spaxel sep- 
arately. The upper panel of Fig. |6] shows the line-of-sight veloc- 
ity (left) and the line FWHM (right) the [Cii] emission line gas 
as inferred with Herschel. The X- and Y-axis are the offsets along 
right ascension and declination relative to the radio core position 
corresponding to the 5 GHz VLBA (very large baseline array) core 



maximum ( [Taylor et al. [2006) 1. These maps have a nominal scale 
of 3" per spaxel. The lower panel of Fig. |6] shows the line-of-sight 
velocity ( left) and the hne FWH M (right) of the HQ'+[Nn] A 6583A 
emission jCanning et al. 



2011 



. The Ha+LNii] X 65 83 A IFU (inte- 
gral field unit) observations were performed with the VIMOS (Vis- 
ible MultiObject Spectrograph) instrument mounted on Melipal, a 
8 m telescope at Paranal Observatory. 

The line-of-sight velocity of the [Cn] emission varies from 
-115 km s"' to 120 km s"'. The gas is blueshifted north of the 
center (represented by the radio core), and redshifted south of 
the center. This is consistent with the velocity map derived from 
Hq'+[Nii] a 6583A emission which shows a similar range in the 
gas velocity and shows that there is gas receding in the south 
and approaching in the north. These results are also in agreement 
with the velocity maps presented by Farage et al.| p010| l using the 
Wide Field Spectrograph (WiFeS) instrument that operates on the 
ANU (Australian National University) 2.3 m telescope at Siding 
Spring Observatory. Fig. |2] shows that the spatial morphology of 
the [Cn] emission closely follows that of the X-ray and optical spi- 
ral filaments. Additionally, the velocity maps show very similar gas 
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kinematics in the infrared and the optical. The velocity and emis- 
sion maps from [Cii] and Ho- suggest a geometry with gas infalling 
on a spiral trajectory starting east of the galaxy center. The trajec- 
tory loops around the galaxy in an anti-clockwise direction, such 
that the infalling gas in the south appears to be receding in the rest 
frame of the galaxy; hence redshifted. The trajectory continues and 
loops around the galaxy towards the center, such that the gas in the 
north appears to be approaching in the rest frame of the galaxy; 
hence blueshifted. 

The velocity dispersion of the [Cn] gas ranges from 
200 km s"' to 600 km s"'. This is slightly larger than the veloc- 
ity dispersion in the H(!'+[Nii] A 6583A gas from the VLT/VIMOS 
observations [50 km s"' to 400 km s"'] and also from the WiFeS 
observations [60 km s"' to 525 km s"']. The larger values for the 
FIR [Cii] data than the optical Ha data is very likely due to beam 
smearing, since the resolution of the Herschel spectral data (11") 
is a lot coarser than that of the optical data ( 0.67"/pix for the VLT 
and 0.5"/pix for WiFeS) 

The morphological and kinematical correlation between the 
far-infrared forbidden line coolant, [Cii], and the optical line fila- 
ments is a key result of this work. This correlation has a profound 
implication, namely that the optical hydrogen recombination line, 
Htt ( Fabian et al.|1982| [Crawford et al.|2005[ [Farage et al.|2010| 



Table 4. Observational constraints. The first column is the quantity, the sec- 
ond column is the measured flux and the third column is the ratio of the 
measured flux to the [Cn] 157.74 /jm flux. 



[Canning et al.||201 

(Dopita et al.||2010[ [Farage et al.||2010[ 

i jCrawford et al.[[2005^ and the far-infrared [Cn] line all have the 
same energy source. 



the optical forbidden lines, [Nn] X 6583A 
the soft X-ray filaments 



5 MODELING THE INTERSTELLAR MEDIUM OF 
NGC 4696 

To understand the complete picture giving rise to emission in 
NGC 4696 at different wavebands ranging from X-ray to radio, we 
performed a detailed investigation of the physical parameters of the 
interstellar medium of the BCG. 

The close spatial correspondence between the [Cn] maps, 
along with the NW-SE elongated dust emission, and the optical Ho- 
maps strongly favours a common heating mechanism for the gas. 
The excitation source that leads to the optical line-emission and 
far-infrared coolants remains to be investigated. The Herschel ob- 
servations presented in this paper shed light on a key ingredient of 
this developing picture. The positive detection of the most luminous 
far-infrared coolant of the ISM, [Cii], along with [Oi] and [Nii], dic- 
tates the presence of photodissociation regions. In this section we 
study the source(s) responsible for heating the photodissociation 
regions (PDRs). 

PDRs form a sharp interface with adjoining Hii/Hi and Hn re- 
gions. Photons from nearby young stars (O- and B-type) with en- 
ergies > 13.6 eV ionize the hydrogen in the molecular gas giving 
rise to an Hii region. This is followed by a PDR where FUV pho- 
tons with energies 6 eV < < 13.6 eV are absorbed by dust grains 
which re-radiate energy in the form of FIR continuum. About 0. 1 % 
to I % of the incident energy of the FUV photons is converted 
into photoelectrons with energies < 1 eV which heat the gas via 
elastic collisions. With CLOUDY we attempted to model an inte- 
grated Hii/PDR cloud, where the output emission originating from 
the cloud also contains contribution from layers beneath the PDR 
surface. 



Quantity 


Value (10"'^ erg s"' cm"^) 


F('i)/^'Cii 


[Cii]" 157 //m 


174.7±3.1 


1 


^FIR 


(1.105±0.119)xlO'' 


63.3;'7,3 


[Oi]'-' 63 urn 


57.6±7.7 


0.33!«;°^ 


[Nuf 122 tim 


24.9±1.7 


13+"'" 


Ha'' 


79.0±0.9 




[Sii]-^ 


<4 


< 0.023 


[On:]" 


<3 


< 0.017 


[OlB]*^ 


< 2 


< 0.011 



" Integrated FIR flux in the ran ge (40-500) / im 
Extinction corrected Ho' from 



Farage et al. 



2010 



This work. [Sii], [Oni] and [OibJ are non-detections. 

5.1 Setting up the CLOUDY simulations 

A detailed modeling of the various observed continuum and line 
emission in NGC 4696 was conducted using the spectral synthesis 
code for modeling the photodissociation regions, CLOUDY (ver- 
sion 08.00 described by Ferla nd et al.[p^98) . CLOUDY models 
the clouds based on conservation of energy through balancing the 
cooling and heating rates. The simulations are self-consistent and 
include chemistry, radiative-transfer and thermal balance. The ba- 
sic ingredients of each simulation comprised an old stellar pop- 
ulation (OSP) and a young stellar population (YSP), wherein the 
Starburst99 synthetic library ( [Leitherer et al.[ 1999{[Vazquez & LeF] 
[therer[2005) was used to model the spectra of both the stellar popu- 
lations. In addition, a background radiation was included compris- 
ing the UVX cosmic background from radio to X-ray (e.g. |Ikeuchi[ 
[& Ostr iker 1986 , Ved el et al.|I994[ l, which also included the cosmic 
microwave background radiation. 

For the YSP, we used the same spectrum as used to derive 
the FUV star formation rate in Sect. [3.5| where the oldest stel- 
lar age was set to 10' yr and the star formation rate was con- 
sidered to be continuous and equal to 0.1 Mg yr"'. For the OSP, 
the oldest stellar age was set to 10 Gyr and the mass was fixed to 
10" Mq. Among the specified model parameters were the normal- 
izations for the OSP, 7/osp> and the YSP, GO, in units of the Habing[ 
1 11968) field, where one Habing is the interstellar radiation field of 
1.6 X 10"^ erg s"' cm"^ in the energy range 6 eV and 13.6 eV. An- 
other parameter was the total hydrogen density (atomic, molecular, 
ionized and other hydrogen-bearing molecules) in units of cm"'' . 

This study and several other studies of NGC 4696 provide ob- 
servational constraints, which may be exploited to constrain the 
physical properties of the ISM of the galaxy. The data explicitly 
used to determine the most likely scenario are listed in Table|4] The 
errorbars on [Cii]/Lpi[( were estimated based on the minimum and 
maximum of the possible values derived from Lfir calculated with 
the 6 Herschel + MIPS 24 yum data points and LpiR calculated with 
the averaged data points in case of duplicate measurements (see 
Sect. [3.3[ ). The lower errorbar on [Oi]/[Cii] was derived by integrat- 
ing the flux over all the spaxels with a SNR> 3 for the [Cii] line and 
a SNR> 2 for the [Oi] line (because of the lower instrument sen- 
sitivity at the wavelength of the [Oi] line emission and a generally 
weaker emission relative to [Cii]). The upper errorbar on [Oi]/[Cii] 
was derived by assigning the spaxels with no [Oi] detection but 
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Table 5. The basic photodissociation region (PDR) model parameters. The geometry is assumed to be plane-parallel. Columns: (1) parameter, (2) symbol, (3) 
the input range and (4) the most likely parameter values. 



Parameter 


Symbol 


Input Range 


Likely Values 


Total Hydrogen Density (cm^^) 


n 


10 to lO*" 


50 to 100 


FUV Intensity Field (Habing*) 


GO 


1 to 10** 


10 to 80 


Extra Heating (erg s"' cm"^) 


^extra 


10-24 iQ-20 


~ 10-22 


Normalization for the OSP* (10"'* erg s"' cm"^ Hz"') 


»70SP 


2.4 to 156.7 


2.4 


Hydrogen Column Density (cm"^) 


Nh 


10" to 1023 


102' 


Metallicity 


z 


1 -4 


2.5 


Nitrogen abundance (relative to Z) 


Zo(N) 


2-2.5 


2 



" 1 Habing = 1.6 X 10-^ erg s"' cm-2, OSP: Old Stellar Population 



0.02 
0.018 
0.016 
0.014 
0.012 

0.01 
0.008 
0.006 
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Figure 7. The metallicity and penetration depth dependency of the line ratios, [Cnl/LpiR (upper left), [Oi]/[Cii] (upper right), HQ'/[Cn] (lower left) and 
[Nn]/[Cn] (lower right). The different lines coiTespond to the different penetration column densities, Nh = lO" cm-2 (red solid hne), 102" cm-2 (cyan dashed 
line), 102' cm-2 (dark blue dotted line), 1022 cm-2 (magenta long-dashed dotted line) and 102^ cm-2 (orange short-dashed dotted line). The other parameters 
have been fixed to GO = 30 Habing, n = 100 cm^^, H^yx^.^ = 10^22 erg s"' cm"' jhe dust and nitrogen abundances are 1.2 times and 2 times the default ISM 
values. The shaded regions correspond to the upper and lower limit of the observed ratios. 



a positive [Cn] detection a bare minimum flux equal to the sensi- 
tivity of the PACS spectrometer in the third order times the chosen 
SNR, and likewise for the [Cn] line. The errorbars on [Nii]/[Cn] and 
Hq'/[Cii] reflect the statistical uncertainties in the measurements. 

In addition, |Kaneda et al.| ( |2005^ found the IRS (Spitzer In- 



frared Spectrograph) spectrum of NGC 4696 to show only one sig- 
nificant polycyclic aromatic hydrocarbon (PAH) feature at 12.7 yum 
with a flux estimate of 14.9 x 10"'^ erg s"' cm"2 pAH features are 
very prominent in the vicinity of young massive stars and so are 
considered as good tracers of star-forming sites. Absence of any 
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strong PAH features is in concert with the fact that both the FIR- 
and FUV-derived star formation rates are low. Also, even though 
the dust-to-gas ratio in NGC 4696 is higher than the nominal range 
(Sect.[34](, PAHs owing to their small size are susceptible to grain 
destruction via physical sputtering or thermal evaporation ( jPwek &] 
|Arendt|1992| l. While PAHs have an important effect on the chem- 
istry of PDRs fBakes & Tielens |(1998 1, due to the lack of any strong 
PAH features in the IRS spectrum we did not consider PAH grains 
in our simulations. 

The geometry was chosen to be plane-parallel by making the 
inner radius much larger than the thickness of the cloud. The in- 
cident radiation from the different stellar populations was input in 
units of flux (erg s"' cm"-). The simulated output fluxes were de- 
termined relative to the [Cn] flux and compared to the observations. 
As a check on the absolute flux levels, we determined the luminos- 
ity of the simulated [Cii] 157.74 yum line by multiplying the flux 
with the surface area of the emitting cloud and scaling it to the dis- 
tance of the Earth. Although the radius of the cloud can be judged 
from the size of the sphere encompassing the [Cn] emission shown 
in Fig.|6]as R ~ 17" (~ 3.5 kpc) (also see Sect[3TTJ, the [Cii] flux is 
not expected to fill up this whole sphere. The volume filling factor 
can be determined by equating the simulated and observed fluxes. 
Note that the PDR model thus constructed is oversimplified in that 
a single uniform PDR cloud has been assumed. In reality, one ex- 
pects a distribution of PDR clouds with different and inhomoge- 
neous incident and emergent radiations. However, as shown below, 
a simple model such as described here is sufficient to study some 
of the basic physical properties of the ISM in NGC 4696. 

X-ray observations of the Centaurus cluster show supersolar 
abundances in the inner 25 kpc, where the metallicity is ~ 2 ( |Gra-| 
|ham et al. '2006 1. Hence we increased the elemental abundances, Z, 
to two to three times the default ISM values. Furthermore, owing 
to a small gas-to-dust mass ratio derived in Sect. |3.4| we increased 
the dust abundance to 1.2 times the default ISM value leading to a 
commensurate decrease in the gas-to-dust mass ratio from ~ 150 to 
~ 70. Another important parameter is the penetration depth of the 
incident radiation into the PDR cloud. The outer radius of the cloud 
was specified in terms of the total hydrogen column density, Nh, at 
the edge of the cloud (the PDR surface). While the observed ratio 
of [Nil] to [Cii] demands relatively low Nh, the observed ratio of 
Hq- to [Cn] demands high Nh. We investigated a range of Nh and 
Z for the most likely values to explain the observed ratios. The de- 
pendency of the ratios on Nh and Z is shown in Fig.|7] Based on the 
results of this investigation, we fixed Z to 2.5 and Nh to 10^' cm"-. 
The various model parameters and their input values can be found 
in Table[5] 

5.2 Photoionization by stars 

The aim of conducting these simulations was to understand the key 
physical ingredients behind the extended [Cn] and Ha filaments. 
Is a young stellar population important? Can the old and young 
stellar populations alone reproduce the observed ratios listed in Ta- 
ble[4]? Do the observational constraints require an additional source 
of heating? In the AGN-heating regulated feedback framework, the 
AGN outflows return a fraction of the accreted energy back to the 
intracluster gas. The radio source can heat the gas either via en- 
ergetic cosmic ray particles or inflating cavities. The other sce- 
nario entailing additional heating revolves around merger-induced 
shocks in the ISM of the BCG (see Sect.|6). 

The young stellar population is important in terms of provid- 
ing FUV photons capable of ionizing carbon to produce the ob- 
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Figure 8. Simulations without any young stellar population produce low 
levels of [Cn] emission in comparison with the observations. As a re- 
sult the [CnJ/LpiR ratio is several orders of magnitude below the observed 
value. The different curves correspond to the different levels of ?/osp, 
the normalization for the old stellar population, increasing top to bottom 
from 2.4 X 10"'* erg s"' cm"^ Hz"' (yellow short-dashed dotted line) to 
156.7 X 10""' erg s"' cm"^ Hz"' (red solid line). 



served [Cn] emission. Without any young stars, the [Cii]/Lfir ratio 
can not be reproduced. This is shown in Fig. [5] where different 
curves represent different values of t/osp- The old stellar popula- 
tion produces optical and UV photons which heat up the dust but 
does not affect the gas much. Hence increasing //osp has an effect 
of increasing LpiR but not the intensity of the gas line emissions, 
which decreases the [Cii]/Lfir ratio. Though elementary, this sim- 
ple exercise underlines the importance of having both the stellar 
populations. Furthermore, the observed [Cii] to LpiR ratio ranges 
from 0.1 % in normal and starburst galaxies to I % in dwarf ir- 
regular galaxies ( [Kaufman et al.|1999[|Luhman et al.|2003^ . |Edge| 
|et al.|p010b|a[ > studied the [CnJ/LpiR ratio in two cool-core BCGs 
and found it to be equal to ~ 0.4% for Abell 1068 and ~ 1.2% 
for Abell 2597. In comparison, the observed ratio for NGC 4696 is 
~ 1.6 %. Thus the ratio is higher than normal and requires mini- 
mum amount of dust heating. While such a large ratio may be rec- 
onciled with the PDRs in low-metallicity galaxies with low dust-to- 
gas mass ratios, NGC 4696 has a very high dust-to-gas mass ratio. 
In order to reproduce the observed [Cn] to LpiR ratio, we fixed the 
770SP to have a bare minimum value of 2.4x 10"'* erg s"' cm"^ Hz"' . 
A lower value of 770SP does not produce any discernible effects on 
the output ratios. 

Following this, we did simulations based on a model contain- 
ing both the old and young stellar populations, and probed n and 
GO space ranging from (10 to 10*) cm"^^ and (10 to 10*) Habing. 
The left panel of Fig. [9] shows the contours of the simulated ratios 
[Oi]/[Cii] and [Cn]/LpiR corresponding to the lower and upper limit 
given in Table |4] (where reproducible). We find that the [Cii]/LpiR 
ratio (dotted blue line) saturates at an upper limit of about 0.014, 
where as the observed ratio is 0.0147. Also, with this model the 
Hq'/[Cii] ratio is greater than 1 for all values of n and GO, where 
as the observed ratio is close to 0.45. Since neither the observed 
[Cnj/LpiR nor the Hf!'/[Cn] can be reproduced, we conclude the 
model containing only old and young stellar populations as the 
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source of heating the PDRs is insufficient to explain the observa- 
tions of NGC 4696. 



5.3 Extra heating 

The gas in PDRs is mostly heated via the dust through photoelec- 
tric heating. In order to boost the [Cii] flux independent of the FIR 
dust emission, a source of heating is required that may directly in- 
put energy to the gas. Thus we included an extra heating term in 
CLOUDY using the parameter "hextra", which specifies a volume- 
heating rate (erg s"' cm"^). The source of this extra heating is un- 
specified and so may be attributed to either shock heating or cos- 
mic ray heating due to the radio source. We investigated models 
containing the OS? and YSP and additional heating ranging from 
10"^* erg s"' cm"^ to 10"^" erg s"' cm"'. The term "hextra" acts 
to raise the temperature of both the gas and the dust independent 
of each other. The standard behaviour of [Cii] is such that for a 
given GO the [Cii] intensity increases linearly up to the critical den- 
sity (n ~ 3000 cm"' for the [Cii] forbidden line) beyond which 
coUisional de-excitation occurs as often as radiation. Hence the 
line intensity above the critical density is not very sensitive to the 
density and increases with it only very slowly. The extra heating 
acts to boost the [Cn] intensity at low densities all the way up to 
n ~ IC* cm"' after which it decreases due to overly suppressed 
radiative de-excitation. The extra heating also increases the dust 
emission; however the increase is roughly independent of the den- 
sity. The over all effect of the heating term on the [ChJ/Lfir ratio is 
such that the ratio increases for low densities, where the increase in 
[Cn] is larger than the increase in Lfir- 

The middle panel of Fig. |9] shows the contours of [Cii]/Lfir 
(solid red), [Oi]/[Cii] (dotted blue) and Hq-/[Cii] (dashed green) for 
^^extra = 10""^ erg s"' Cm"'. Of the different levels of //extra in- 
vestigated, this was the optimal value leading to most coherent re- 
sults. For example, //exc^ = 10"^' erg s"' cm"' fails to reproduce 
the observed [Oi]/[Cn] ratio and //extra = 10"^' erg s"' cm"' just 
barely touches the observed [CnJ/LpiR ratio. We found that while 
//extra could reproducc the right levels of [Cii]/LpiR, Hcf/[Cn] and 
[Oi]/[Cii], the [Nii]/[Cn] ratio was too low compared to the obser- 
vations. The FIR [Nil] line emission (the first ionization potential of 
nitrogen is 14.5 eV) arises from warm ionized medium only, where 
as [Cn] emission can arise from both neutral (PDRs) and ionized 
media. The fraction of [Cn] emission arising from ionized media 
is highly uncertain, ranging from 25 % up to more than 50 % (e.g. 
|Heiles|1994[|Oberst et al.|2006[[Amnestad & Emery 2003 1. How- 
ever, since our model includes a PDR adjacent to an Hii region, 
the total [Cn] emission includes contributions from both, neutral 
as well as ionized media, and no correction factor is needed to esti- 
mate the fraction of [Cn] emission arising from the PDR alone. The 
right panel of Fig.|9]corresponds to the same setting as for the mid- 
dle panel, only the nitrogen abundance, Zo(N), has been increased 
by a factor of two over the ISM value. The contours corresponding 
to the [Nii]/[Cii] ratio (dotted-dashed mustard) now overlap with 
the contours of the [Cii]/LpiR, [Oi]/[Cii] and Hq'/[Cii] ratios. A ni- 
trogen overabundance from the PDR analysis is consistent with the 
enhancement in nitrogen abundance required based on the optical 
line ratios parage et al.|2010l > and X-ray spectral fitting l |Sanders| 
let al.|2008| >. The fact that the analyses of the NGC 4696 observa- 
tions in X-ray, optical and far infrared wavebands, all require a ni- 
trogen overabundance and by roughly the same factor is reassuring 
and lends support to the above PDR model. 

The emanating [Cn] flux from the cloud is about 
1 X 10"' erg s"' cm-. The size of the cloud that results in the 



observed [Cn] flux is such that, Fq^^ ^j^j, = Fq^-^ mod ^ (''/^l)'^' 
where Fq^^ qJ^^ is the observed [Cii] line flux and Fq^^ mod 
the modelled [Cn] line flux, r is the "effective radius" of the cloud 
and /)l is the luminosity distance to NGC 4696. The radius of the 
cloud thus calculated is r = 0.5 kpc. Considering that the [Cn] 
flux extends over a ~ 3.5 kpc large region implies a volume filling 
factor of ~ 3 X 10"'. 

These PDR simulations help us to narrow down the most prob- 
able heating scenario such that it simultaneously reproduces all the 
observed line and continuum ratios. The most likely parameter val- 
ues are given in the last column of Table |5] The right panel of 
Fig.|9]is based on a model that contains photoionization from old 
and young stellar populations, and an additional source of heating. 
Note that photoionization from young stars is still of paramount 
importance. A model containing only extra heating and no young 
stars produces very little [Cn] because of severe lack on FUV pho- 
tons. The model containing stars and "extra-heating" best explains 
the observations yielding a density, n, in the range of about a few 
tens to hundred per cm' and a radiation field, GO, in the range ten 
Habing to ~ 80 Habing. This range of GO corresponds to a star 
formation rate of (0.02 to 0.13) Mq yr"' for a cloud of an effec- 
tive radius of 0.5 kpc. The SFR thus estimated agrees very well 
with the limits derived from the FIR and NUV/FUV observations 
(Sect.[331l. 



6 DISCUSSION 

The Centaurus cluster of galaxies can be classified as a strong cool- 
core cluster based on its short gas cooling time (< 0.5 Gyr) and 
a peaked X-ray surface brightness profile. The mass condensation 
rate derived from the far-infrared and UV observations is small (~ 
O.lMo yr"') as compared to that predicted based on the cooling- 
flow model (~ 25 Mq yr"'). The X-ray core luminosity at 2.5 %/?5oo 
corresponding to a region of radius 30 kpc is about 1 .3x lO"" erg s"' 
( [Sanders et al.|2008} . This is the amount of power (radiative losses) 
that needs to be fed back to the ICM over that scale for the system to 
remain in equilibrium. One of the potential feedback mechanisms 
is through the expanding radio plasma that has dug out multiple 
cavities in the intracluster medium of NGC 4696 at radii < 6 kpc. 
The power in these radio bubbles assuming a buoyancy timescale 



is estimated to be about 0.7 x 10"" erg s"' ( [Rafferty et al.|2006 i 



Repetitive inflation of such radio bubbles also creates sound waves 
which generate a sound-wave power of about 0.5 x 10"" erg s"' 
(Sanders et al.||2008} . These two mechanisms together suffice to 
balance the X-ray radiative losses and can keep the central region 
of Abell 3526 in quasi-equilibrium. 

The optical and far-infrared line emissions clearly show an in- 
falling gas flow of spiral pattern. There are two possibilities for the 
origin of the gas and dust. The first possibility is that this material 
represents the cooling intracluster medium flowing toward the cen- 
ter of the dark matter potential well (which in the case of Centau- 
rus coincides with the BCG). The fact that the soft X-ray emission 
(0.5-1) keV traces the spiral structure seen in Ha and [Cn] does 
indeed suggest a link to the gas cooling out of the ICM. The rela- 
tive increase in the X-ray surface brightness along the filaments is a 
natural outcome of this scenario. jFabian et aLHTOTl) examined the 
Chandra X-ray emission from the Northern filament in the Perseus 
cluster and provided another outlook on the excitation mechanism 
of the filaments. Their results indicate that some of the soft X-ray 
emission may be due to charge exchange as the hot ionized gas pen- 
etrates and mixes with cold gas. Although the cooling flow model 
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Figure 9. Simulations containing botli old and young stellar populations. The dust abundance is set to 1.2 times the ISM value and the metallicity is set to 
2.5 times the ISM value. Shown are the contours corresponding to the lower and upper observed limits (where reproducible) for [Cii]/Lfir (solid red line), 
[Oi]/[Cn] (dotted blue Une), Hq'/[Cii] (dashed green line) and [Nn]/[Cn], (dashed-dotted mustard line). Left: Heating by stars only. Middle: Heating by stars 
and an additional source, such as cosmic-rays, shocks etc. Right: Nitrogen overabundance of a factor of two. 



does not dictate a coherent directional gas flow to the center, such 
as that seen in NGC 4696, the spiral morphology may be a natural 
consequence of the global velocity field of the gas surrounding in 
NGC 4696. This may be the first direct mapping of a cooling flow 
in a galaxy cluster. 

The second possibility is that the filaments and the dust lane 
in NGC 4696 have a galaxy infall origin. This is based on a merger 
hypothesis proposed by [Sparks et al.| ^1989^ and more recently 
by Farag e~et al.| p010| l. Under this hypothesis, the BCG is un- 
dergoing a minor merger with an infalling neighbouring dust-rich 
galaxy. The presence of two nuclei in HST imaging of the center 
of NGC 4696 (j Laine et al.|2003) lends some support to the merger 
hypothesis, although the nature of the second nucleus is not well 
established. A merger hypothesis also explains the prominent dust 
lane seen with the HST. While [Sparks et al.| jl989| l favour heat 
transfer via conduction from the hot X-ray emitting ICM to the 
cool dust-rich infalling gas, |Farage et al.| ( (2010[ > favour shock heat- 
ing generated through the dissipation of the orbital energy of the 
infalling cloud. Assuming a shock luminosity of 2 x lO'*^ erg s"' 
( [Parage et al.|201()) , the shock volume luminosity over a region of 
radius 0.5 kpc is 1.3x 10"^^ erg s"' cm"^. This is in good agreement 
with the level of the extra heating input to the simulations above. 

Irrespective of the origin of the cold gas and dust, the detailed 
PDR modeling of the observed quantities arising from this flow 
(described in Sect. |5j shows photoionization from stars is an im- 
portant requisite. However, there is an additional source of energy 
input required to reproduce the observables. This additional source 
of heating could either represent the energy deposited by the AGN 
through the methods described above or shock heating due to a 
minor-merger. The additional heating amounts to < 10 % of the X- 
ray core luminosity and could be another contributing factor to the 
total heating budget. 

The non-radiative energy sources, in principle, could be in- 
put into the simulations either in the form of extra heating (shock 
heating or reconnection diffusion, see [Fabian et al.[2011[ for latter), 
such as that considered in these CLOUDY simulations, or cosmic 
rays. [Ferland et al.[ l [2008) and [Ferland et al.[ ( [2009 1) have shown 
that while it is difficult to dilTerentiate between the two processes, 
certain optical lines could help discriminate between the two. The 
extra heating case in those studies has primarily been associated 
with magneto-hydrodynamic (MHD) waves such as those observed 
in the ISM of our galaxy. The dissipation of MHD wave energy can 



heat the gas; however, the heating simply adds to the thermal en- 
ergy of the gas so that the velocity distribution remains Maxwellian. 
On the other hand, the cosmic ray heating by high-energy particles 
can both heat and ionize the gas. Both the processes heat the gas 
but the cosmic ray (ionizing-particle) case produces a population 
of first and second ions by non-thermal coUisional ionization. In 
the extra-heating case these ions only occur when the gas is warm 
enough for collisional ionization equilibrium to occur. This is the 
main distinction between the two cases, extra heating and cosmic 
ray heating. 

,Ferland et al.j ( [2009^ studied the Horseshoe region of the 
Perseus galaxy cluster and used the infrared and optical line inten- 
sities to distinguish between extra heating and cosmic ray heating. 
They found that both the heating cases match the observations to 
within a factor of two for the majority of the lines. There are a few 
discriminant lines, such as the optical emissions [Hei] A 5876 A 
and [Ne m] A 3869 A and the infrared emission [Ne ii] A 12.81 yum 
which show a few orders of magnitude difference and indicate that 
the cosmic ray may be a better agent for heating and ionizing the 
gas (also see Donahue et al. 2011 1. The mid-infrared Spitzer IRS 
spectrum of the off-center regions of NGC 4696, in fact, contains 
the [Ne ii] A 12.81 yum line ( [Johnstone et al.|2007| l. This along with 
the detections of ro-vibrational H2 lines in off-nuclear regions in 
NGC 4696 and NGC 1275 provided some of the motivation for 
the heating models investigated by [Ferland et al. (20091. However, 
previous studies did not consider the lines observed with Herschel, 
which provide a crucially important perspective on the physical 
model of the ISM. 

Though our model contains additional energy explicitly in the 
form of extra heating, we speculate the ionizing-particle model 
would also have fitted the available data, as has been shown for 
Perseus. In order to discern the two heating sources, deep optical 
observations are needed to complement the Herschel data. 



7 SUMMARY 

We made far-infrared (FIR) observations of the brightest galaxy, 
NGC 4696, of the Centaurus galaxy cluster with the Herschel tele- 
scope to better understand the cooling and heating of the intraclus- 
ter medium. 

• We have detected FIR coolants in NGC 4696, which include 
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extended [Cii] 157.74 /jm, marginally extended [Oi] 63.18 yum and 
unresolved [Nii] 121.90 yum line emissions. 

• We have detected far-infrared dust emission from the BCG at 
70 yum, 100 yum, 160 /jm, 250 //m, 350 yum and 500 fim. A spec- 
tral energy distribution (SED) fitting of the dust emission reveals 
a high-mass cold component around 19 K and a low-mass warm 
component around 60 K. 

• Using the dust mass calculated from the SED fitting and the 
non- detection of CO ( |0'Deaet al.[T994) , we derived an upper limit 
on the gas-to-dust mass ratio of 125. This makes NGC 4696 a 
galaxy with one of the lowest gas-to-dust mass ratios. 

• The star formation rate derived from the integrated FIR lumi- 
nosity is about 0.13 Mq yr"'. We derived similar upper limits from 
HST far-ultraviolet and GALEX near-ultraviolet observations. 

• The extended [Cn] emitting gas shows remarkably similar spa- 
tial morphology and kinematics as the optical Ha emitting and X- 
ray gas. This implies a common heating mechanism of the gas. 

• We envisage the FIR and optical emissions as arising from 
photo-dissociation regions (PDRs) adjacent to ionized regions. 
From a detailed modeling of such an integrated PDR, we conclude 
that in addition to heating via stellar photoionization, an additional 
non-radiative heating is required. The most likely model yields a 
total hydrogen density in the range (50 to 100) cm"' and a FUV 
intensity field in the range (10 to 80) Habing. 
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